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ABSTRACT 

The ejecta of the Cas A supernova remnant has a complex morphology, consisting 
of dense fast-moving line emitting knots and diffuse X-ray emitting regions that have 
encountered the reverse shock, as well as more slowly expanding, unshocked regions of 
the ejecta. Using the Spitzer 5-35 /an IRS data cube, and Herschel 70, 100, and 160 /an 
PACS data, we decompose the infrared emission from the remnant into distinct spectral 
components associated with the different regions of the ejecta. Such decomposition 
allows the association of different dust species with ejecta layers that underwent distinct 
nuclear burning histories, and determination of the dust heating mechanisms. Our 
decomposition identified three characteristic dust spectra. The first, most luminous one, 
exhibits strong emission features at ~ 9 and 21 /an, and a weaker 12 /an feature, and 
is closely associated with the ejecta knots that have strong [Ar II] 6.99 /an and [Ar III] 
8.99 /an emission lines. The dust features can be reproduced by magnesium silicate 
grains with relatively low MgO-to-Si02 ratios. A second, very different dust spectrum 
that has no indication of any silicate features, is best fit by AI2O3 dust and is found in 
association with ejecta having strong [Ne II] 12.8 /an and [Ne III] 15.6 /an emission lines. 
A third characteristic dust spectrum shows features that best matched by magnesium 
silicates with relatively high MgO-to-SiCU ratio. This dust is primarily associated with 
the X-ray emitting shocked ejecta and the shocked interstellar/circumstellar material. 
All three spectral components include an additional featureless cold dust component 
of unknown composition. Colder dust of indeterminate composition is associated with 
[Si II] 34.8 fim emission from the interior of the SNR, where the reverse shock has 
not yet swept up and heated the ejecta. The dust mass giving rise to the warm dust 
component is about ~ 0.1 M 0 . However, most of the dust mass is associated with the 
unidentified cold dust component. Rs mass could be anywhere between 0.1 and 1 M 0 , 
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and is primarily limited by the mass of refractory elements in the ejecta. Given the 
large uncertainty in the dust mass, the question of whether supernovae can produce 
enough dust to account for ISM dust masses in the local and high-z universe remains 
largely unresolved. 

Subject headings: dust, extinction — infrared: ISM — ISM: individual (Cassiopeia A) 
— supernova remnants 


1. INTRODUCTION 

Interstellar dust models that fit the average interstellar extinction curve, the diffuse infrared 
emission and scattering, and abundances constraints (Li &; Draine 2001; Zubko, Dwek, & Arendt 
2004; Brandt & Draine 2012) employ a very limited variety of dust compositions, consisting primar- 
ily of polycyclic aromatic hydrocarbons (PAHs), graphite or amorphous carbon, and astronomical 
silicates. Yet observations of specific interstellar dust sources, meteoritic inclusions, and interplan- 
etary dust particles (IDPs) reveal a significantly richer variety of dust compositions. Examples are 
magnesium sulfides (MgS) found in the atmospheres of AGB stars near the Galactic Center (Chan 
et al. 1997), and presolar grains of SN or stellar origin, such as silicate carbide (SiC), silicon nitride 
(Si 3 N 4 ), and aluminum, calcium- and titanium-oxides (e.g. Zinner 2008). Their absence in inter- 
stellar dust models suggests that they are not required for fitting the various manifestation of dust 
in the general ISM, either because of their low abundance relative to silicates and carbonaceous 
dust, or due to the fact that they may have been processed in the ISM. 

Supernovae can be important sources of interstellar dust. They produce all the refractory 
elements needed for the formation of dust, and their ejecta largely retain the compositional in- 
homeogeneity of the progenitor star. They may therefore be sources of dust with unusual chemical 
compositions. Furthermore, SN are the drivers of the chemical evolution in galaxies, and therefore 
potentially the most important sources of interstellar dust. In fact, SN are the only source of newly 
condensed dust in high redshift galaxies, before low mass stars evolved off the main sequence. De- 
termining the mass and composition of SN condensed dust is therefore key for understanding the 
origin, evolution, and processing of dust in galaxies. 

The Cas A remnant is an ideal object for studying the composition and abundance of dust 
that formed in the ejecta of core collapse SNe. The SN was not definitively recorded at the time 
that it occurred. Studies of the expansion of the SNR estimate that the explosion was in the year 
1681Y19 (Fesen et al. 2006). Yet fortunately, light echoes of the Cas A SN have allowed studies 
of this old event with modern instruments. Such observations have revealed that Cas A was a 
Type lib SN (Krause et al. 2008) with distinct asymmetry in its explosion (Rest et al. 2011). 
Dynamical and compositional asymmetries are still imprinted on the Cas A SNR today, but the 
dominant structure of the Cas A SNR is characterized by a clear distinction between the forward 
shock sweeping up the interstellar (and/or circumstellar medium) and the reverse shock through 
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which the SN ejecta is expanding. 

The ejecta consists of three main components: the first, containing most of the mass, consists 
of a low density phase that is heated by the reverse shock to X-ray emitting temperatures (> 10 6 
K). The second component consists of dense fast-moving knots (FMKs) that have gone through 
the reverse shock, and are radiatively cooling by line emission at optical and IR wavelengths. A 
third component comprises ejecta that has not yet encountered the reverse shock, and is primarily 
heated by the ambient radiation within the SNR interior. 

In this paper we revisit the analysis of the mid-IR spectra of the dust in Cas A. The infrared 
emission arises from dust residing in the circumstellar/interstellar (CS/IS) medium that is shocked 
by the advancing SN blast wave, from dust in the diffuse X-ray emitting and dense line emitting 
knots that encountered the reverse shock, and from dust in yet unshocked regions of the ejecta. 
Our main goal is to decompose the IR emission into its different spatial components in order to 
determine the mass, composition, and the specific regions of the ejecta in which the dust was 
formed. Our approach will therefore provide important information on the physical processes that 
facilitate or inhibit the nucleation of dust in the different layers of the ejecta, the resulting dust 
composition, and the dust heating mechanisms that give rise to the IR emission. Our approach is 
significantly different from previous ones (Ennis et al. 2006; Rho et al. 2008) which only grouped 
the IR emission into distinct spectral components, with no relation to the nature of the ejecta from 
which they originated. 

To decompose the IR spectrum into the different components of the ejecta, we use multi 
wavelength observations of the remnant as proxies for their composition and spatial distribution. 
We use the radio synchrotron emission to represent the shocked CS/IS medium, the Chandra X-rays 
emission to represent the shocked diffuse ejecta, and the Spitzer fine-structure IR emission lines to 
represent the dense FMK. The description of the data sets used in the analysis and the construction 
of the different spectral templates representing the different components of the ejecta are described 
in §2. In §3 we describe the decomposition of the dust into the different spectral templates, and 
derive the IR spectra associated with each template. In §4 we derive the minimal dust mass and 
composition required to fit the IR emission associated with each template. The discussion in §5 
assembles the results to create a clearer picture of just what can and cannot be constrained from 
the analysis of the mid-IR emission of Cas A. The overall results of the paper are summarized in 
§ 6 . 


2. DEFINING THE SPECTRAL/SPATIAL TEMPLATES OF THE EJECTA 

2.1. Data Preparation 


Work began with the low resolution IRS spectral data cubes used by Rho et al. (2008). 
There are 4 cubes, generated from the 1st and 2nd order spectra in short and long wavelength low 
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resolution modules. At each wavelength, the data were convolved to the resolution at the longest 
IRS wavelength, 38.33 //rn, using kernels derived from Tiny Tim/ Spitzer (STINYTIM) 7 PSFs 
according to Gordon et al. (2008). The three shorter wavelength cubes were then reprojected to the 
same scale and orientation as the longest wavelength (SL1) cube. The data are then combined into a 
single IRS low-resolution cube covering 5-38 /im with a fixed spatial resolution. Weighted averages 
were used at the wavelengths that overlap between the different spectral orders and modules. 

For the purpose of analyzing the dust emission, we created a continuum spectral cube by 
replacing all the emission lines with a smooth polynomial fit to the continuum on either side 
of each line. This procedure is performed independently for each spatial pixel of the cube and 
each line. It avoids any need to fit the lines themselves, which can be found at various velocity 
components at each pixel because of the high expansion velocities of Cas A. The lines that were 
thus removed are listed in the row and column headings of Table 1. Finally, to remove remaining 
artifacts and improve the signal to noise, the continuum cube was smoothed in wavelength to a 
spectral resolution of R = A/AA ~ 20, and trimmed to a maximum wavelength of 35 //rn. 

Longer wavelength information was obtained by using Herschel PACS observations at 70, 100, 
and 160 //rn, with the appropriate spatial convolution of the 70 and 100 /im images whose native 
resolution is better than the IRS at 38 //rn. At 160 //m the SNR is substantially confused by ISM 
emission. This confusion was reduced by subtracting an empirically scaled version of the SPIRE 
250 //rn image in which the synchrotron emission of the SNR was subtracted by extrapolation from 
the 350 and 500 gm bands. The subtraction is an improvement, but is still imperfect, in part due 
to the lower spatial resolution of the longer wavelength SPIRE data. 

The synchrotron emission is removed from the final Spitzer+ Herschel continuum cube by 
subtraction of the VLA radio template (below) extrapolated assuming a power law spectrum S u ~ 


2.2. Spatial Templates 

A set of templates are needed to identify the spatial distribution of various physical components 
of the ejecta. The primary source for templates are the line emission maps that are generated by 
subtracting the continuum cube from the full data cube, and then integrating over the [-6000,-1-8000] 
km s _1 velocity range for each of the emission lines in the IRS spectrum. Most of these lines are 
associated with the fast-moving knots of ejecta. These knots tend to be O-rich, but there are some 
variations in composition that tend to be correlated with spatial and kinematic differences. Because 
of the large velocity range in Cas A, we did not disentangle the closely spaced [O IV] 25.89 //m and 
[Fe II] 25.99 //m lines. 


'http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/contributed/general/stinytim/ 
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Several other spatial templates were also tested. These are intended to trace the presence of 
dust in a wider range of environments than the relatively cool and dense line-emitting knots. (1) 
A radio template was derived from 6 cm VLA image of Cas A (DeLaney 2004). The region of the 
forward shock is traced more clearly in the radio than at most other wavelengths, and the bright 
radio knots in the reverse shock region are a different population than line-emitting ejecta knots. 
(2) and (3) X-ray Fe and Si templates were derived from Chandra observations (Hwang et al. ????), 
integrated over the widths of the 6.4 keV Fe and 1.7 keV Si K a lines. The X-ray continuum is 
included in these integrations. These templates trace ejecta in regions of lower density and much 
higher temperature. This include the bulk of the mass of the ejecta of the SN. (4) the IRAC 4.5 
/mi image lies outside the range of the IRS spectral coverage and contains both a line emission 
component and a synchrotron component. The synchrotron emission was subtracted using a scaled 
version of the IRAC 3.6 fi m image which is dominated by synchrotron emission and has only weak 
line emission. This also over-subtracts stellar sources in the image. Thus regions where negative 
values resulted were replaced by nearby background values. (5) The IRAC 8 /.mi image is dominated 
by the [Ar II] 6.99 /mi emission, but could potentially be a more sensitive tracer of this emission 
than the IRS line maps. The background ISM emission at 8 /.mi was removed by subtracting a 
scaled version of the 5.8 /mi IRAC image, where the SNR emission is relatively weaker. As with the 

4.5 /.mi template, the over-subtracted stellar sources were replaced by the local background levels. 
Each of these templates is convolved and reprojected to match the IRS continuum data cube. 

This full set of templates would be expected to contain some degeneracies, where different lines 
(especially from the same species or element) would be tracing the same physical component of 
the SNR. To select a useful subset of the possible templates, we calculated the linear correlation 
coefficient for each pair of templates. These are listed in Table 1. The table shows which templates 
are well-correlated with others. Based on these correlations we sorted the templates in to similar 
groups and selected one template to represent each group. Images of the selected templates are 
shown in Figure 1. 

There are several points of interest in the correlations listed in Table 1. Even though the [Ar II] 
and [Ar III] emission are very well correlated (as expected), we chose to keep both these templates 
in order to monitor for changes in the dust that might be related to the ionization state of the gas. 
It was not expected that highly ionized species such as [S IV] would correlate so well with [Ar II], 
This correlation may become significantly weaker if higher spatial resolution data were available. 
Correlations involving the [Fe II] line are weaker than most due to the low signal to noise ratio 
of this line. The [Ne V] line is also relatively weak, yet surprisingly it correlates better with the 
[Ar II] group than the [Ne II] and [Ne III] lines. Despite an absence of Ne lines in the band, the 

4.5 /m i IRAC /im emission is somewhat better correlated with the Ne group templates than the 
Ar group. In particular, the “Ne crescents” noted by Ennis et al. (2006) are present in the IRAC 

4.5 /tm image. Although the [S III] 33.48 /an correlates moderately well with the [S III] 18.71 /an 
line, it does not correlate as well with the other templates of Ar group. This is an indication of 
variation in the S III line ratios, which would indicate large scale variations in the density of the 
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ejecta. The [O IV] + [Fe II] template is distinct from the others, but we did not use it because as a 
blend of two lines it is likely to represent a mixture of different dust types and environments. The 
X-ray Si template is in fact distinct from the X-ray Fe template, but it was not used because initial 
analysis using only the X-ray Fe template did not indicate the presence of any residual emission 
that would have been compensated for by the addition of an X-ray Si template. 


3. DERIVING THE INFRARED SPECTRA OF THE DIFFERENT TEMPLATES 

3.1. Characteristic Template Spectra 

One of the primary goals of this analysis was to find the continuum spectra that are associated 
with zones of various physical parameters and compositions throughout the SNR. To identify the 
spectra associated with each of the six emission templates shown in Fig. 1, we superimposed 
all the templates and for each one we identified the zones where its relative emission was the 
dominant component. Because of the strong and widespread [Ar II] emission, in some case the 
zones necessarily included a significant [Ar II] emission as well. The zones selected for each of the 
templates are outlined in Figure 2. 

For the selected pixels, i , of each zone, we performed a linear fit to the spectra, D(i, A): 

6 

D(i,X) = Y / S j (X)T j (i) + C( A) (1) 

3 = 1 

where Tj(i ) are the 6 selected emission templates. The derived parameters are (7(A), a constant 
term (which should be small) to account for errors in the background levels, and Sj( A), the spectra 
of each of the emission components. The characteristic spectra derived are shown in Figure 3. 

When these derived spectra are multiplied by the spatial templates over the entire SNR (not 
just the selected zones where each was derived) and the resulting model is subtracted from the 
data cube, we found that there was one region of particularly strong residual emission, consistently 
positive at all wavelengths. This region is in the south central portion of the SNR. One additional 
zone was created to cover this region (the “South Spot” in Fig. 2), and the spectrum extracted 
from this zone is shown in Fig. 3. 


3.2. Spatial Distribution of the Characteristic Spectra 

The final step in the decomposition of the continuum spectral cube is to repeat the sort fitting 
given by Equation (1), but this time we use the seven Sk( A) spectra as known quantities, and solve 
for the spatial distributions of each component over the entire SNR. 

7 

D(i,A) = ]Ts fc (A)T fc (i) + C(i) 

k= 1 


(2) 
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In this fitting, a constant term C(i) is still present to account for background errors, but now 
the constant is a function of position rather than wavelength. The spatial templates derived here, 
Tfc(i), should be similar to the initial templates Tj(i), but would also identify regions where dust 
with a similar emission spectrum is found, but is superimposed on other emission components. The 
derived spatial distributions of the different characteristic spectra are shown in Figure 4. 


3.3. Residuals 

The total emission of each of the seven components is calculated by integrating the spectral 
templates, Sj( A) (Fig. 3), over the regions of the corresponding spatial templates, X/ (Fig. 4). 
Figure 5 shows the seven Cas A spectra with their actual intensities (not normalized). This com- 
parison shows that in the ejecta, the Ar II dust component is dominant at wavelengths from 8 to 
35 /mi. The dust associated with the radio emission is assumed to be swept up ISM dust, and is 
generally < 1/3 of the luminosity of the ejecta dust. The sum of the all the spectra is also shown 
as a “fit” to the total Cas A spectrum derived without any decomposition into components. The 
general agreement at 10 - 35 /im indicates that there are no major missing components. At shorter 
and longer wavelengths, the comparison is not as good. This is primarily due to difficulties in the 
background subtraction when integrating over the entire area of the SNR. 

Images (at selected wavelengths) of the differences between the total Cas A emission and the 
sum of the components are shown in Figure 6. Examination of these images and the full set of 
residuals at all wavelengths shows that at any chosen locations the residuals can be both positive 
and negative as a function of wavelength. Prior to the inclusion of the South Spot component, there 
was a conspicuous positive residual at nearly all wavelengths at that location. The variations of the 
residuals indicate that locally there are more subtle variations of the emission in the SNR than can 
be traced by our derived spectral templates. Such variations would be expected from small changes 
in the dust temperature that would be expected if the physical conditions that prevail where each 
type of dust is found are not strictly uniform across the SNR. 


4. DUST COMPOSITIONS FOR DIFFERENT TEMPLATES 

Six of the seven characteristic dust spectra show strong emission at wavelengths from 10 - 
40 |jm. The Spitzer IRS spectra across these wavelengths reveal distinctive features in each of 
these spectra. The features generally appear to be suggestive of silicates with peaks of varying 
sharpness at ~ 10 and ~ 20 /im. In order to identify the particular type of silicates and/or other 
grain compositions that produce each of the characteristic spectra, we assembled a set of grain 
emissivities from published optical constants (Table 2), and fitted each template’s spectrum, Sk{ A), 
as the weighted, M n> j, sum of different grain emissivities, /c n (A), applied to blackbody emission at 
different temperatures, B(Ti ): 


- 8 - 
2000K 

S k ( A) = E E M nji K n (\)B(Ti) . (3) 

n Ti= 20K 

The MPFIT code (Markwardt 2008) was used to determine coefficients M n i needed to fit each 
of the normalized spectrum with the constraint that M n ^ > 0. For the purpose of fitting, the 
spectrum uncertainties were taken to be <75 = 0.25 with a floor of as > 0.015. In the fitting, the 
dust temperature is constrained to lie on a grid of values ranging logarithmically from T = 20 to 
T = 2000 with Alog 10 (T) = 0.05. Only components with intensities exceeding 10 -4 times the peak 
of the fitted spectrum are kept in the final result. 

For each spectrum, an initial set of fits was performed using only one dust composition, n, 
at a time. No single composition provided a good fit to the spectrum. The fits were repeated 
using all possible combinations of pairs of compositions (i.e. the left summation in Equation 3 
is over n = 0, 1). Some combinations explored in this way result in substantial improvements in 
the fits. Usually the improved fits involved the pairing of a dust composition with strong spectral 
features with another composition that is relatively “featureless.” Thus, we also performed a set 
fits using 3 compositions, but with the constraint that one was required to have a “featureless” 
emission spectrum. Since all the featureless compositions produce similar results, we chose the “ac” 
amorphous carbon composition for the fixed third component in these tests. 

The quality of the one and two component fits as measured by x 2 is shown in Figures 7-12. 
The left panel of each figure shows a matrix of the x 2 values, with generalized axes. This is useful 
to get an overview of which pairs of components yield provide the best fits. The general groupings 
are based on chemical and/or spectral similarity. The right panel of each of the figures displays the 
same x 2 values but now plotted as a function of the second composition and color coded according 
to the first composition. This plot is more useful for illustrating the quantitative changes in y 2 and 
allows specific components to be identified. The results for single composition fits correspond to 
the upper envelope of the curves in these plots (or the diagonal of the matrix in the left panels). 
Examples of the best fitting 2- and 3-composition models, as well as contrasting models that also 
provide reasonably good fits, are shown in Figures 14-34. The results for the 7 different spectra are 
discussed below. 


4.1. Ar II dust 

The emission associated with Ar II exhibits the sharpest peaks at 9 and 21 /in 1 and seems 
to have a weaker peak at ~12 /tm. The astronomical silicate of the typical ISM does not have 
silicate peaks sharp enough to produce the observed features. Magnesium silicates (particularly 
Mgo. 7 Si0 2 .7, characterized by MgO / Si0 2 = 0.7) have peaks that are in roughly the right locations, 
but their shapes are not a good match to the data. However, when two-composition fits are 
performed, the addition of a second component lacking strong spectral features (e.g. graphite, 
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amorphous C, or AI 2 O 3 ), can smooth soften and smooth the silicate peaks to provide a much 
better fit to the spectrum. 

The best 2-component fit (Figure 14) employs Mgo. 7 Si 02.7 to fit the 9 and 21 /zm silicate 
peaks, and graphite as the featureless component that weakens the relative strength of the peaks. 
The fit to these peaks is very good, although the weaker 12 //rn peak is not fit at all. Despite the 
overall good fit, the temperature components needed for the Mg silicate are somewhat suspicious. 
Components at 71 — 63K are needed to fit the 21 //rn feature, but are too cold to contribute to the 
9 /zm feature. Conversely a 500 K component produces the 9 /zm peak while contributing little 
to the 21 /zm peak. Small, stochastically heated grains might be expected to be heated to this 
range of temperatures, but it is problematic that there is no indication of grains at intermediate 
temperatures (70K< T < 400K). This suggests that the emissivities used here are not a true 
physical representation of the Ar II dust. The actual dust must either have a relatively stronger 9 
/zm peak so that it can produce the observed spectrum from cooler dust, or it must have a somewhat 
lower emissivity at roughly 10-20 /zm so that the presence of grains at intermediate temperatures 
would be required. 

Figure 15 shows an alternate 2-conrponent fit where the same Mg silicate is paired with a 
nonstoichiometric spinel [Mg^Al^CH; y = (8 — 2x)/3]. These spinels (especially with a low Mg/Al 
ratio) have a spectral peak that improves the fit to the observed 12 /zm feature, at the expense of 
poorer fits to the other silicate peaks. The maximum temperature required for the silicate grains 
is much lower here, but the temperature of the spinel reaches 2000K, which is the maximum value 
allowed. This high temperature is driven by two influences. The first is to add the 12 /zm feature 
to the spectrum without making significant contributions at longer wavelengths. The second is to 
provide the (uncertain) continuum level at < 8 /zm. The fits are not seriously worse if the maximum 
temperature is restricted to a lower value (e.g. 400K). 

As suggested by the previous 2 -component fits, the best 3-component fit combines the Mgo.ySi 02.7 
with a featureless amorphous C (to weaken the silicate peaks) and a nonstoichiometric spinel to 
provide the 12 /zm peak (Fig. 16). However there is still a tendency in this fit to use distinctly 
different temperature and/or composition components to fit each different spectral feature. The 
spinel material that produces the 12 /zm peak is still required to be at extremely high temperatures. 
Alternate sources of the 12 /zm feature are SiCD and SiC. For both of these materials, the spectral 
features are too sharp and slightly too blue to provide a good match when the mass absorption 
coefficients, k(A) are calculated from Mie theory for spherical grains. However, calculations assum- 
ing a continuous distribution of ellipsoids (CDE) for particle shapes, result in spectral features that 
are broader and redder. The dust temperatures required to produce the features are also much 
more reasonable. Figures 17 and 18 show the fits when CDE calculations are used for a third 
composition of SiC >2 or a-SiC. Si 02 provides some contribution to all three of the observed peaks. 
SiC only provides the 12 /zm peak. 
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4.2. Ar III dust 

The Ar III dust spectrum is very similar to that of the Ar II dust. It was derived separately 
to check if there might be an identifiable shift in dust temperature or composition in these closely 
related regions. Not surprisingly, the combinations that produced good fits to the Ar II spectrum 
also produce good fits to the Ar III spectrum (Figures 19-16). This increased strength of the 9 /rrn 
peak relative to the 21 /im peak would usually suggest warmer dust temperatures in association 
with the Ar III. However, because the dust features arise largely from separate components of the 
models, the temperature of the components are largely unchanged, but the relative mass in the 
hotter components is higher for the Ar III than for the Ar II spectrum. 

For the Ar III spectrum we consistently find that the amorphous C (ac) provides a slightly 
better fit than the graphite that was somewhat preferred for Ar II spectrum. The Ar III spectrum is 
also slightly better fit when hibonite rather than the spinels is used as a third component, although 
if CDE calculations are applied, Si02 and SiC can again provide equally good fits as the third 
component. 


4.3. Ne II dust 

The dust associated with Ne II exhibits an unusually smooth spectrum. This spectrum excludes 
significant amounts of warm silicates, but can be fit by compositions with relatively featureless 
spectra, e.g. graphite, and amorphous carbon. The best fits are found by combining featureless 
dust with AI2O3. The latter component has a very broad emission peak between 10 and 20 //rri that 
improves the fit. Figure 22 shows the fit for an amorphous carbon (be) - AI2O3 mixture. The fit 
suggests that some fraction of the AI2O3 component is fairly warm (T ~ 150 K), while the brightest 
carbon dust has a temperature of ~ 75 K. A significant cool (~ 40 K) dust component helps to fit 
the Herschel data at 70 - 160 /iin. The model attributes this to AI2O3, but this identification is 
degenerate because neither of these compositions have any characteristic spectral features at these 
wavelengths. The composition of very cool components cannot be identified very conclusively, 
because the shape of the spectra are hardly affected by composition. 

Figure 23 shows the best fit model that does not contain AI2O3. In this case a nonstoichiometric 
spinel provides a relatively small perturbation to an otherwise smooth amorphous C spectrum. 
Amorphous C picks up the intermediate temperature (~ 135K) components that account for the 
bulk of the 10 - 20 //m emission. 

Figure 24 shows the best 3-conrponent model for the Ne II dust. The addition of Ti02 nomi- 
nally improves the fit over the best 2-conrponent model, but the changes to the spectrum are too 
minor to suggest that any third component is warranted. 
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4.4. X-Ray Fe dust 

The spectrum of dust associated with the X-ray Fe emission does exhibit apparent silicate 
peaks at ~ 10 and 20 /an. The 10 //m peak is relatively weak, suggesting a lack of hot silicate 
grains. As with the Ar II dust, mixing silicate and more featureless emission (e.g. amorphous 
C) improves the quality of the fit. The preferred silicates here tend to have higher Mg0/Si02 
ratios than those that fit the Ar II and Ar III dust. Figure 25 shows the best 2-composition fit. 
The temperatures for both components are warm, T ~ 100 — 112 K, with additional hot and cold 
components to account for the short and long wavelength emission. Figure 26 shows the best 
2-component fit that is obtained without the use of silicate compositions. This marginal fit uses 
Mgo.6Feo.4O to adjust the shape of the model spectrum near 20 /im, but contains no replacement 
for a 10 /mi silicate feature. The best 3-conrponent fit (Fig. 27) merely adds PAH + emission to 
reproduce a bump at 8 /mi and sharpen the spectral feature at ~ 10 /im. The overall change in y 2 
is small. 


4.5. South Spot dust 

The South Spot spectrum looks rather different than the X-ray Fe spectrum, yet it tends to be 
best fit by similar dust compositions. Here the temperatures of the components that produce the < 
20 /mi emission are higher than those for the X-ray Fe dust. Figure 28 shows the best 2-component 
fit which uses similar compositions as the X-ray Fe spectrum (cf. Fig. 25). Figure 29 shows a 
marginal fit obtained without the use of silicates, and thus without any component that provides a 
10 /mi feature in the spectrum. The best 3-component fit (Fig. 30) uses nonstoichiometric spinel 
to make minor adjustments to the shape of a good 2-component model. 


4.6. Radio (ISM) dust 

As for the X-ray Fe and South Spot spectra, the spectrum associated with the radio emission 
tends to be best fit with a mix of silicate and featureless dust. Figure 31 shows the best 2- 
composition fit. Figure 32 shows a marginal 2-composition fit, in which Mgo.1Feo.9S serves as the 
“featureless” component across the 10-20 /mi part of the spectrum. The best 3-conrponent fit (Fig. 
33) is not significantly different from the best 2-component fit, using two featureless components 
instead of one. Figures 34 shows the best results when dust compositions are restricted to the 
astronomical silicate, graphite, PAH and PAH + , which are commonly used to fit general interstellar 
material (e.g. Zubko et al. 2004). The shape of the astronomical silicate’s 10 /mi emission peak is 
not as sharp as the observed spectrum. Cool PAH + emission provides sufficient adjustment to the 
shape of the silicate spectrum such that the model does not require and graphite or neutral PAH 
components. 
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4.7. Si II dust 

The Si II dust spectrum essentially contains only the Herschel PACS measurements at 70, 100, 
and 160 /in i, with only upper limits on the emission at shorter wavelengths. Because of this lack 
of detailed spectral information, the spectrum was only fit with single composition dust models. 
The Si II spectrum has a relatively sharp peak compared to a black body spectrum. The best fit 
is provided by CaAl^Oig which has a broad emission feature at ~ 80 /im making it an especially 
good fit to the data. Figures 35 and 36 show the best fit spectrum (CaAl^Oig), and a fit using 
a more typical Mg silicate (Mg 2 SiC> 4 ). Compositions that have a steeper spectral index at long 
wavelengths provide better fits, but the relatively good fit for Ti02 (3) (rutile) is an artifact of 
an unphysically steep extrapolation (A -3 5 ) implied by the measured optical constants, limited to 
< 35 fun. For the Si II spectrum, any models that are within a factor of 4 of the minimum y 2 are 
deemed acceptable. This includes most dust compositions, and only excludes compositions that 
have strong features near ~ 40 /tm that would have exceeded the upper limits provided by the 
Spitzer IRS data. 


5. DISCUSSION 

The identifications of possible dust compositions in the different environments in Cas A are 
summarized in Table 3. A single dust composition can never provide a good fit to the observed 
spectra, except for the Si II spectrum which has only upper limits at < 70 /tin. In general, 2 
compositions are sufficient to get acceptable fits to the spectra. The Ar II and Ar III spectra are 
the only ones that show significant (though small) benefit for the addition of a third composition. 
Among the 7 characteristic dust spectra examined, the results can be grouped into 3 different 
families of dust. 

The first family is found in association with the Ar II and Ar III emitting ejecta, which are 
distributed widely across the SNR. This family consists of the “lightest” Mg silicate (i.e. lowest 
MgO/Si 02 ratio) in combination with one or more other compositions that have a featureless 
spectrum. The Mgo.7Si02.7 is a good fit to the 9 and 21 /an peaks in the observed spectrum, 
but only if a featureless dust composition is also present to reduce the apparent strength of these 
features. However, the Ar dust spectra also contain a weaker 12 pm feature which is not accounted 
for by any Mg silicate. Nonstoichiometric spinel with low Mg/Al ratios can provide a feature at 
approximately the correct wavelength, but only when the dust grains are extremely hot, such that 
longer wavelength spinel features are relatively faint. A better explanation for the 12 pm feature 
may be provided by Si02- The spectrum of SiC >2 exhibits exhibits all three peaks seen in Ar II 
and Ar III spectra, although they are are significantly sharper and slightly bluer than the observed 
ones, and even in combination with other materials the fits are not very good. However, this is 
when Mie theory is used to calculate the absorption cross sections assuming small spherical grains. 
Using a continuous distribution of ellipsoids (CDE) approximation instead, broadens and shifts the 
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Si 02 features to be a better match as a third component. Jager et al (2003) point out that the 
12 /mi SiC >2 feature disappears in Mg silicates when Mg 0 /Si 02 > 0.5. Therefore, it seems likely 
that the dust associated with the Ar emission is a Mg silicate - silica mix with Mg 0 /Si 02 < 0.5 in 
combination with a separate featureless dust component. SiC can also provide the 12 /t m feature, 
but again only if CDE calculations are applied to this component. If SiC is present, this would be 
the only direct evidence of carbon-bearing dust. Table 4 lists the possible origins for the 12 /mi 
feature of the Ar dust. 

The strong 21 /mi peak in this family has been the hallmark of Cas A IR spectra since it was 
first observed using KAO and ISO. On the basis of those data, the peak was suggested to arise 
from Mg protosilicate (Arendt et al 1999). Subsequent analysis by Douvion et al. (2001) modeled 
an ISO spectrum as MgSiOs, Si 02 and AI 2 O 3 with the weak 12 /mi feature largely produced by the 
AI 2 O 3 as proposed for the feature in “Spectrum 2” of Douvion et al. (1999). Ennis et al. (2006) 
using Spitzer IRS data noted the distinction of several different dust spectra in different parts of 
Cas A and referred to this as the “Strong 21 /mi” spectrum, and Rho et al. (2008; calling it “21 
/mi peak dust”) modeled it as Mg protosilicate and MgSiC >3 with secondary components of SiC> 2 , 
FeO, FeS, Si, AI 2 O 3 or FeO, SiC> 2 , FeO, FeS, Si, Fe. 

The second dust family is associated with the Ne II emission, which is especially prominent in 
two opposing “Ne crescents” (Ennis et al. 2006; Smith et al. 2009) in the N and S parts of the 
SNR. These morphological features are evident in the distribution of the Ne II dust shown in Figure 
4. This dust has a very smooth spectrum that does not suggest any silicate material. The best 
fits to the spectrum are found with AI 2 O 3 in combination with other featureless dust. Featureless 
dust alone can provide moderately good fits, but a broad asymmetric feature in the 10 - 20 /mi 
portion of the AI 2 O 3 absorption cross section seems to match the observed spectrum especially well. 
Alternately, nonstoichiometric spinel can provide the needed emission at 10 - 20 /mi, although the 
spinel emissivity has more detailed substructure that is not evident in the observed spectrum. 

This second family corresponds to the “weak 21 /mi” components noted by Ennis et al. (2006) 
and Rho et al. (2008) which they also associate with relatively strong Ne emission. However, the 
fact that they see even a weak 21 /mi peak in their spectrum suggests that it is a mixture of what 
we identify as very distinct Ar and Ne dust families. As in our fitting, Rho et al. (2008) fit this 
spectrum with hot and cool featureless dust (C glass) and with an intermediate temperature AI 2 O 3 
components. They included other components to add the weak 21 /mi peak that appears in their 
spectrum. Douvion et al. (1999) had also noted an anti-correlation between the 9 /mi silicate 
emission and Ne II and Ne III emission. 

The third dust family is associated with the X-ray Fe emission and the South Spot. This family 
also seems to match the dust associated with the radio emission, which is largely dust that has 
been swept up from the interstellar or circumstellar medium by the forward shock. The primary 
component in this family is one of the “heavier” Mg silicates (i.e. with high MgO/Si 02 ratios) 
or MgFe silicate: Mg 2 SiC> 4 , Mg 2 . 4 SiC> 4 . 4 , MgFeSiC> 4 . Other Mg silicates with MgO/Si 02 > 1 can 
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provide acceptable fits, but the silicates with lower ratios that were needed for the Ar dust are not 
suitable here because of the changing placement and shape of the 9 and 21 features. 

This dust family is matched by the “Broad” component identified by Ennis et al. (2006) and 
the “Featureless” component modeled by Rho et al. (2008). In both cases this component is not 
associated with Ar or Ne emission lines, just as we also find. Rho et al. (2008) modeled this 
component as MgSiCU, FeS, and Si combined with AI2O3, IVfeSiOj and/or Fe. Using 5-17 gm ISO 
data, Dovioun et al. (1999) extracted a spectrum (“Spectrum 3”) from a region that should match 
our Radio spectrum. Despite the limited wavelength coverage, they also found that the spectrum 
could be fit with Draine & Lee astronomical silicate (but without graphite) at T ~ 105 K, as we 
confirm in Fig. 34. 

The Si II dust has no significant emission across the 5-40 /um wavelength range of the IRS. The 
three HerschelPACS measurements at 70, 100, and 160 are insufficient to constrain the composition 
of the associated dust. The Si II dust may belong to one of the above families, or it may be an 
entirely different composition. 

The modeled dust temperatures and compositions allow the determination of the dust mass 
in each of the components. Generally the dust mass is dominated by a warm component (60 
- 130K) that also produces the bulk of the luminosity. However the X-ray Fe spectrum is an 
exception, because the Herschel PACS data at 70-160 //rri are elevated relative to the shorter 
wavelength emission, which requires an additional cold component that likely dominates the mass. 
The composition is uncertain because emissivities are smooth and the spectral resolution is poor 
at the long wavelengths. This is uncertainty is worse for the Si II dust which is only detected 
at > 70 /un. Its temperature is relatively well constrained, but its composition, and therefore its 
mass, is not. If the dust is composed of Mg silicates, then the total mass is < 0.1 M 0 , but this 
value can be either lower or much higher if other compositions are appropriate. The derived dust 
masses averaged over all models that fit within a factor of 2 of the minimum y 2 and are dominated 
by Mg silicates (or AI2O3 for the Ne II spectrum) are listed in Table 3. Figures 37-43 plot the 
total dust mass for all 2-component models as a function of y 2 with color coding to indicate 
compositions and temperatures. The “good” models are within a factor of 2 of the minimum \' 2 
(to the left of the dashed line). The figures show that although there are correlations between 
composition and temperatures, the derived mass is more strongly dependent on the composition 
than the temperature. 


6. SUMMARY 

We have decomposed the IR spectrum from Cas A into contribution from distinct regions of 
the remnant, comprising: the shocked CS/IS medium into which the remnant is expanding, the hot 
X-ray emitting gas and fine-structure IR line emitting regions that encountered the reverse shock, 
and unshocked regions of the ejecta. We then calculated the mass and composition of the dust 
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associated with each region of the ejecta. The methodology and results of our paper can be briefly 
summarized as follows: 

• We first identified a set of spatial templates representing the IR emitting regions of Cas A. 
These were defined by the fine structure IR lines from Ar II, Ar III, Ne II, and Si II, by the 
X-ray (Fe line), and by the radio synchrotron emission. The spatial templates used for the 
decomposition are shown in Figure 1. 

• We then identified spatially distinct zones, shown in Figure 2, in which a given template emis- 
sion was dominant. We then extracted the dust spectra from these zones using the procedure 
described in Section 3. The resulting IR spectra were assumed to represent the different 
spatial templates, and are presented in Figure 3. Figure 4 shows the spatial distribution of 
the dust associated with each of the spectral templates. 

• We compiled an extensive list of dust compositions with known optical properties to calculate 
their contribution to the IR emission from each spatial templates. The dust species fall into 
8 broad categories: silicates, protosilicates, silica, silicon carbide, carbon and metallic iron 
dust, aluminum oxides, oxides, and sulfides. The 59 dust compositions used in the analysis 
are listed in Table 2 

• For each region we fitted the spectrum with all possible combinations of 2 distinct dust species, 
equaling a total of 59x58/2 combinations. Figures 7-13 depicts the 59x59 matrices of the \' 2 
values from all the possible combinations. Figures 14-36 show how the addition of additional 
dust species affect the value of % 2 for various fits to the template spectra. 

• The primary dust component could be readily identified in regions that had strong IR dust 
features. The Ar II and Ar III ones exhibited strong IR features at 9 and 21 /mi. Magnesium 
silicates, characterized by a Mg0/Si02 stochiometric ratio of 0.7 (i.e. Mgo.7Si02.7) provided 
the best fit to the spectra of these regions. The Ne II region was best fit with AI2O3 dust, 
and the X-ray Fe region, which exhibit peaks at ~ 10 and 20 /mi was best fit with silicates 
having a higher Mg0/Si02 ratio of 2.4 (i.e. Mg2. 48104.4). A region called the South Spot had 
similar composition. The dust composition from the radio synchrotron emitting region was 
consistent with that expected from interstellar dust. The IR spectrum from the Si II emitting 
region arises from a cold dust component that is only seen at long wavelengths. Since very 
few dust compositions have distinguishing features at these wavelengths, the composition of 
this cold dust could not be determined. 

• The secondary dust component needed to improve the fit for most regions is not uniquely 
identified. Table 3 lists the primary dust composition in each region in bold, and all possible 
secondary dust species that they need to be paired with. The secondary dust species generally 
have a featureless dust spectrum. 

• Altogether, the IR spectrum of each of the seven regions in Cas A can be fit with two dust 
species. Some species are common to more than one region, so a minimal number of only 
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4 dust species: Mgo. 7 Si 02 . 7 , Mg 2 . 4 SiO 4 . 4 j AI 2 O 3 , and amorphous carbon, are required to fit 
the entire spectrum of Cas A. These compositions suggest that the seed dust particles that 
formed the more complex species are MgO, Si 02 , AI 2 O 3 , and carbon. 

• The mass of dust in the different zones is presented in Figures 37-43. The mass of dust is 
about 0.04 Mq in the shocked CSM/ISM and ejecta regions, and <0.1 M 0 in the unshocked 
ejecta characterized by the Si II emission. This dust mass is similar to that derived by Rho 
et al. (2008) from the Spitzer data, and the ~ 0.08 M 0 derived for the cold unshocked ejecta 
derived from the Herschel and Akari data (Barlow et al. 2010, Sibthorpe et al. 2010). 

• The total mass of newly-synthesized dust (~ 0.1 — 0.14 M 0 ) is not sufficient to account for 
the large amount of dust detected in high-redshift galaxies, unless the dust is not destroyed 
by the expanding SN blast wave in the ISM (Dwek, Galliano, &; Jones 2007). 


This work is based on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
Support for this work was provided by NASA Program NNH09ZDA001N-ADP-0032. This research 
made use of Tiny Tim/Spitzer, developed by John Krist for the Spitzer Science Center. The Center 
is managed by the California Institute of Technology under a contract with NASA This research 
has made use of NASA’s Astrophysics Data System Bibliographic Services. 
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Fig. 1. — Spatial templates used for the spectral decomposition. The [Ar II] 6.99 /mi. [Ne II] 12.81 
/jin, [Si II] 34.8 n m, [Ar III] 8.99 fim lines are derived from the IRS data cube. The X-ray Fe 
emission is from Hwang et al. (????) and the VLA 6 cm radio map is from DeLaney (2004). 
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Fig. 2.— Depiction of the zones for which the different spectra were extracted. These zones corre- 
spond to the regions where the templates shown in Figure 1 are dominant (apart from the nearly 
ubiquitous Ar II). The additional “South Spot” zone was identified as the primary region where 
there is relatively strong continuum emission that is not traced by any of the spatial templates. 
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Fig. 3.— The continuum spectra extracted for each of the spectral zones defined in Figure 2. The 
solid line is the IRS data. The three points at 70, 100, and 160 /im are broadband Herschel PACS 
data. 
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Fig. 4. — Distribution of the dust associated with each of the extracted spectral types shown in 
Figure 3. 
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Fig. 5. — (left) Comparison of the relative intensities of each of the spectral types with the total 
emission of Cas A. The gray line labeled “Fit” is the sum of all the component types, (right) 
Comparison of the total Cas A emission with the sum of all the ejecta dust components and the 
ISM dust component. The ISM dust component is that associated with the radio emission of the 
forward shock. 


Declination Declination Declination 


-24- 



23 " 23 '" 45 ° 30 ° 15 ' 

Right Ascension 


00 ° 23 23 45 ° 30 ° 15 ° 

Right Ascension 


Right Ascension 


Fig. 6. — Images of the residual emission at select wavelengths in the IRS spectral cube after 
subtraction of the sum of the dust components shown in Figures 3 and 4. 
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Fig. 7 . — (left) A map of x 2 for 2 -composition models of the Ar II dust spectrum. Darker shading 
indicates lower values. The map shows that the lowest values of x 2 are achieved when one of the 
silicates (second row from bottom) is paired with the relatively featureless dust such as C, Fe, AI2O3. 
Good fits are also obtained by pairing this silicate with nonstoichiometric spinels (MgAl204, i.e. 
oxides) with low Mg/Al ratios, (right) In this alternate depiction of x 2 for 2 -composition models, 
each line plotted shows x 2 for a given “primary” dust component (as coded by line/label colors) 
when paired with each possible “secondary” component (as listed in order along the abscissa). 
Thus each plotted line corresponds to a row in the matrix shown to the left. Here one can see more 
specifically that 40 < x 2 < 80 for Mgo.ySi02.7 when paired with the relatively featureless graphite, 
amorphous C, Fe, AI2O3, Fe3C>4, and FeS, or with the nonstoichiometric spinels with low Mg/Al 
ratios. 
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Ar III 



Fig. 8. — x 2 as in Fig- 7, but for the Ar III dust spectrum. Because the Ar III dust spectrum is 
very similar to the Ar II dust spectrum, y 2 for each pair of compositions is very similar in both 
cases. 



Fig. 9. — x 2 as in Fig. 7, but for the Ne II dust spectrum. The best models contain either porous 
or compact AI2O3 paired with any of the more featureless dust components. 
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Fig. 10. — x 2 as i n Fig- 7, but for the X-ray Fe dust spectrum. This spectrum is best fit by silicates 
with high Mg 0 /Si 02 ratios and paired with the relatively featureless compositions. 



Fig. 11. — x 2 as i n Fig. 7, but for the South Spot dust spectrum. The results here are very similar 
to those for the X-ray Fe spectrum, but with somewhat less contrasts between the best fits and the 
more marginal ones. 
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Radio 



Fig. 12. — x 2 as i n Fig. 7, but for the Radio dust spectrum. The best fits here are for silicates 
with high Mg 0 /Si 02 ratios paired with the more featureless compositions. 



Fig. 13. — X' 2 as in Fig. 7, but for the Si II dust spectrum. Only single composition fits were 
performed because of the lack of detailed spectral constraints. The best fits are for compositions 
that have steeper long wavelength emissivities, and for hibonite (CaAl^Oig) which has a broad 
feature at 80 gm. 
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Total Model 
71 K Mg 07 SiO 27 
1 78K graphite 
63K Mg 07 SiO 27 
502K Mg 07 SiO 27 
1 1 2K graphite 



Fig. 14. — The best 2-composition fit to the Ar II spectrum includes Mgo.7Si02.7 and a featureless 
dust composition. 
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63K Mg 07 SiO 27 
224K Mg 07 SiO 27 
71 K Mg 053 AI 231 0 4 
200K Mg^SiO^ 
2000K Mg 07 SiO 27 
126K Mg 0 53 AI 231 O 4 
2000K Mg 053 AI 231 0 4 
80K Mg 0 53 AI 2 31 0 4 
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Fig. 15. — An alternate 2-composition fit to the Ar II spectrum. Here the nonstoichiometric spinel 
provides the unusual bump at 12 /jrn. 
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Ar II 

Total Model 

71 K Mg 07 SiO 27 

142K graphite 

1 59K graphite 
399K Mg 07 SiO 27 
200K graphite 
71 K Mg 0 73 AI 217 0 4 
356K Mg 07 SiO 27 
2000K Mq n 7 ,AI, 17 0 4 
2000K graphite 



Fig. 16. — The best 3-component fit to the Ar II spectrum is a good match to all of the features 
at 9, 12 and 21 /jrn. However, the 12 /jrn feature is provided by implausibly hot nonstoichiometric 
spinel. 
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Ar II 

Total Model 

71 K Mg 07 SiO 27 

1 59K Si0 2 (CDE) 

200K graphite 
1 1 2K graphite 
71 K Si0 2 (CDE) 
1 78K Si0 2 (CDE) 
399K Mg 07 SiO 27 
224K graphite 
356K Mg 07 SiO 27 



Fig. 17. — A very good 3-conrponent fit to the Ar II spectrum can also be obtained with SiC >2 
as the third component if its mass absorption coefficient is calculated via CDE rather than Mie 
theory. SiC >2 has features at each of the observed peaks. 
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Total Model 
71 K Mg 07 SiO 27 
159K graphite 
142K graphite 
89K a SiC 0 (CDE) 
448K Mg a7 Si0 27 
399K Mg 07 SiO 27 
200K graphite 
80K a SiC 0 (CDE) 
224K graphite 
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Wavelength (^m) 


Fig. 18. — A very good 3-component fit to the Ar II spectrum can also be obtained with SiC as the 
third component if its mass absorption coefficient is calculated via CDE rather than Mie theory. 
SiC contributes only at 12 /mi, but may also provide a partial match of the 21 /mi feature (Jiang 
et al. 2005). 
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Fig. 19. — The best 2-composition fit to the Ar III spectrum makes use of the 12 /mi feature of 
nonstoichiometric spinel. 
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Ar III 

Total Model 
71 K Mg 07 SiO 27 
159K amorph. C (ac) 
50K Mg 07 SiO 27 
63K Mg 07 SiO 27 
564K Mg 07 SiO 27 
502K Mg 07 SiO 27 
45K Mg 07 SiO 27 
142K amorph. C (ac) 



Fig. 20.— Alternate 2-component fits to the Ar III spectrum using featureless compositions fit 
nearly as well, but do not match the 12 /i m peak. 
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Ar III 

Total Model 
71 K Mg 07 SiO 27 
63K Mg 07 SiO 27 
1 59K amorph. C (ac) 
142K amorph. C (ac) 
448K Mg a7 Si0 27 
80K CaAl 12 0 19 
399K Mg 07 SiO 27 
200K CaAl 12 0 19 
224K CaAl 12 0 19 
45K CaAl 12 0 19 
40K CaAl 12 0 19 
2000K amorph. C (ac) 



Fig. 21.— The best 3-component fit to the Ar III spectrum uses CaAl^Oig (hibonite) to fit the 
12 /jrn peak. 
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80K amorph. C (be) 
159K Al 2 0 3 (porous) 

1 42K Al 2 03 (porous) 

7 1 K amorph. C (be) 
31 7K amorph. C (be) 
45K Al 2 0 3 (porous) 
40K Al 2 0 3 (porous) 



Fig. 22.- The best 2-component fits to the Ne II spectrum use AI 2 O 3 in combination with a more 
featureless component. No silicates are a good match for this very smooth spectrum. 
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Total Model 
1 26K amorph. C (be) 
142K amorph. C (be) 
63K amorph. C (be) 
7 1 K amorph. C (be) 
142K Mg 0 53 AI 2 31 0 4 
356K amorph. C (be) 
399K amorph. C (be) 
25K Mg 0 53 AI 2 3 ! 0 4 



Fig. 23. — An alternate 2-component fit to the Ne II spectrum can be found using nonstoichiometric 
spinel instead of AI 2 O 3 . 
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Fig. 24.- The best 3-component fit to the Ne II spectrum adds Ti02, but the improvement in the 
fit is very small. 
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Fig. 25.— The best 2-component fit to the X-ray Fe spectrum uses a Mg silicate in combination 
with a featureless composition. 
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Fig. 26.— This alternate 2-component fit (without silicates) to the X-ray Fe spectrum is about a 
factor of 2 worse than the best fit, and is thus only marginally acceptable. 
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Fig. 27. — The best 3-component fit to the X-ray Fe spectrum uses ionized PAHs to match the 
small (and possibly spurious) bump at 8 ^ m. 
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Fig. 28. — The best 2-component fits to the South Spot spectrum are similar to those for the X-ray 
Fe spectrum, but involved higher dust temperatures for the dominant components. 
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Fig. 29. — An alternate 2-component fit to the South Spot spectrum again shows that reasonable 
(though worse) fits can also be obtained without the use of silicates. 
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Fig. 30. — The best 3-component fit to the South Spot spectrum uses nonstoichiometric spinel as 
a third component to make relatively minor adjustments to the model. 
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Fig. 31. — The best 2-component fit to the Radio spectrum uses graphite and a Mg silicate. 
Combinations of other silicates with other featureless components are nearly as good. 
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Fig. 32. — This alternate fit to the Radio spectrum shows a marginal model using a sulfide com- 
position instead of a one of the more featureless dust compositions. 
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Fig. 33.- The best 3-component fit to the Radio spectrum is only very slightly better than the 
2-conrponent fit. 
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Fig. 34.- The fit of the Radio spectrum using the standard astronomical silicate, graphite, PAH, 
and PAH + combination. Only two of the four components are actually needed for the fit since the 
PAH + component is sufficient to provide a relatively featureless component at > 20 /tm and the 
PAH emission band at ~ 8 /tm. 


Intensity (normalized) 


-50- 



Si II 

Total Model 
32K CaAl 12 0 19 
28K CaAl 12 0 19 
25K CaAl 12 0 19 



Fig. 35. — The best 1-component fit to the Si II spectrum uses CaAl^Oig, which happens to have 
a broad peak in its emissivity at ~ 80 /jm 
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Fig. 36. — This alternate fit to the Si II spectrum uses Mg 2 SiC >4 which has a more typical power 
law emissivity with an index of —2 at long wavelengths. 
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Fig. 37. — Total dust masses for Ar II models. In the top panel the points are color coded by the 
composition that provides the largest mass component. In the bottom panel the masses are color 
coded by the temperature of the dominant mass component. Only models to the left of the dashed 
line (within a factor of 2 in x 2 °f the best model) are deemed acceptable. These good models 
indicate that ~ 10~ 2 M 0 of dust, mostly silicates at ~ 70K is associated with the Ar II emission 
in Cas A. 
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Fig. 38.- Total dust masses associated with Ar III emission. Results are similar to the models for 
Ar II, but the total mass is only half as large. 
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Fig. 39.— Total dust masses associated with the Ne II emission. Here the mass is dominated by 
one of the more featureless compositions or AI 2 O 3 , and the derived total mass depends on the 
long wavelength emissivity for the dominant component. The good models with total masses > 1 
Mq are dominated by Fe dust which has a very low long wavelength emissivity, and thus requires 
extreme dust masses to produce the observed emission. 
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Fig. 40.— Total dust masses associated with X-ray Fe emission. Among the good models, the 
range in mass is caused by the variation in long wavelength emissivity between different dominant 
compositions. 
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Fig. 41.— Total dust masses associated with the South Spot emission. As indicated in Fig. 11, 
there are many different possible compositions that can yield good fits to this spectrum. Thus 
there is a wide spread in mass, depending on what composition is dominant. 
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Fig. 42. — Total dust mass associated with the Radio emission. Again for the best models, there is 
a relatively broad range in mass that is dependent on the composition of the dominant component. 
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Fig. 43.— Total dust mass associated with the Si II emission. Here models within a factor of 4 
of the minimum \ 2 are deemed acceptable, rather than the factor of 2 used in other cases. This 
is because the limited data appear to have underestimated uncertainties and have a coincidentally 
good fit to one dust composition (CaAl^Oig) that has a sharp feature at ~80 /xrn. While all 
compositions must be cold to fit the data, the composition itself, and thus the mass, is not well 
constrained. 
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Note. — Spatial templates in boldface were selected as the basis for further analysis. 


Table 2. Dust Emissivities 
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Composition 
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Note 

Reference 
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Draine & Lee (1984) 
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Mgo. 7 Si0 2 .7 
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Jager et al. (2003) 
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MgSi0 3 

0.2-10000 


Jager et al. (2003) 
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MgSi0 3 

0.1-100000 
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MgSi0 3 
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Jager et al. (1994); Dorschner et al. (1995) 
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Jager et al. (2003) 
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Mg 2 Si04 
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Jager et al. (2003) 

7 
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MgFeSi0 4 
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Jager et al. (1994); Dorschner et al. (1995) 
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Mg 2 . 4 Si0 4 .4 

0.2-8200 


Jager et al. (2003) 

Protosilicates 
10 Ca Protosilicate 

7.6-25 
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Dorschner et al. (1980) 

11 

Ca Protosilicate 

7.7-25 

B - 450C 

Dorschner et al. (1980) 

12 

Ca Protosilicate 
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C - 695C 

Dorschner et al. (1980) 

13 
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Dorschner et al. (1980) 

14 

Fe Protosilicate 

8-40 

E - 490C 

Dorschner et al. (1980) 

15 

Mg Protosilicate 

8-40 

F - 485C 

Dorschner et al. (1980) 

Silica 

16 

Si0 2 

0.1-500 

am 

Kozasa 

17 

Si0 2 

5-500 


Philipp (1985) 

18 

SiO* 

5.6-520 

extrap, at A > 65 pin 

Rinehart et al. (2011) 

SiC, Si 
19 a SiC 

0.1-1000 

21 pm : Qabs/a = 0 

Jiang et al. (2005) 

20 

a SiC 

0.1-1000 

21 pm : Qabs/a = 100 

Jiang et al. (2005) 

21 

a SiC 

0.1-1000 

21 pm : Qabs/a = 10 3 

Jiang et al. (2005) 

22 

a SiC 

0.1-1000 

21 pm : Qabs/a = 10 4 

Jiang et al. (2005) 

23 

P SiC 

0.1-1000 

21 fim . Qabsj O' — 0 

Jiang et al. (2005) 

24 

P SiC 

0.1-1000 

21 pm : Qabs/a = 100 

Jiang et al. (2005) 

25 

P SiC 

0.1-1000 

21 f-LIH . Q abs 1 & — 10 

Jiang et al. (2005) 

26 

P SiC 

0.1-1000 

21 pm : Qabs/a = 10 4 

Jiang et al. (2005) 

27 

Si 

0.02-247.7 

am 

Kozasa 

28 

Si 

0.07-433.3 

cr 

Kozasa 

C, Fe 

29 

Meteoritic Diamond 

0.02-110 


Braatz et al. (2000); Mutschke et al. (2004) 

30 

PAH 

0.01 -9400 


Dwek 

31 

PAH+ 

0.01 9400 


Dwek 

32 

Graphite 

0.01 -9400 


Draine & Lee (1984) 


Table 2 — Continued 


n 

Composition 

A Range (/z m) 

Note 

Reference 

33 

Amorphous Carbon 

0.01-9400 

ac 

Dwek 

34 

Amorphous Carbon 

0.01-9400 

be 

Dwek 

35 

Carbon 

0.1-3000 

glass 

Kozasa 

36 

Fe 

0 . 1-100000 


Kozasa 

A1 Oxides 




37 

A1 2 0 3 

7.8-500 

porous 

Begemann et al. (1997) 

38 

ai 2 o 3 

7.8-200 

compact 

Begemann et al. (1997) 

Oxides 




39 

MgAl 2 04 

2-10000 

natural 

Fabian et al. (2001) 

40 

MgAl 2 O 4 

2-6800 

natural, annealed 

Fabian et al. (2001) 

41 

Mg 1 . 01 Al 1 . 99 O 4 

1.67-6825 

nonstoich. spinel 

Zeidler et al. (2011) 

42 

Mg 0 . 94 Ab. 04 O 4 

1.67-6825 

nonstoich. spinel 

Zeidler et al. (2011) 

43 

Mgo.82Al 2 .i 2 04 

1.67-6825 

nonstoich. spinel 

Zeidler et al. (2011) 

44 

Mgo. 73 Ab. 17 O 4 

1.67-6825 

nonstoich. spinel 

Zeidler et al. (2011) 

45 

Mgo. 53 Ab. 31 O 4 

1.67-6825 

nonstoich. spinel 

Zeidler et al. (2011) 

46 

Mgo. 43 Ab. 38 O 4 

2-10000 

nonstoich. spinel 

Zeidler et al. (2011) 

47 

CaAli20i9 

2-10000 

hibonite 

Mutschke et al. (2002) 

48 

Ti0 2 (1) 

2-5843 

anatase 

Posch et al. (2003), Zeidler et al. (2011) 

49 

Ti0 2 ( 2 ) 

2-5843 

brookite 

Posch et al. (2003), Zeidler et al. (2011) 

50 

Ti0 2 (3) 

0.47-36 

rutile 

Posch et al. (2003), Zeidler et al. (2011) 

51 

CaTi0 4 

2-5820 

perovskite 

Posch et al. (2003), Zeidler et al. (2011) 

52 

Mgo. 6 Feo. 4 O 

0.2-500 


Henning et al.(1995) 

53 

Mgo. 1 Feo. 9 O 

0.2-500 


Henning et al.(1995) 

54 

FeO 

0.2-500 


Henning et al.(1995) 

55 

Fe 3 0 4 

0 . 1-100000 


Kozasa 

Sulfides 




56 

Mgo. 9 Feo. 1 S 

10-500 


Begemann et al.(1994) 

57 

Mgo. 1 Feo. 9 S 

10-500 


Begemann et al.(1994) 

58 

FeS 

0 . 1-100000 


Kozasa 

59 

SiS 2 

13-60 


Begemann et al.(1996) 


Table 3 . Best 2 -Composition Fits of Cas A Dust 


Spectrum 

Compositions 1 

Total Dust Mass (M 0 ) 2 

Ar II 

Mg 0 . 7 SiO 2 .7 + graphite, ac, FeS, C glass, Fe, be, AI2O3, nonstoichiomentric spinel, Fe30 4 

0.01 ±0.002 

Ar ill 

Mg 0 . 7 SiO 2 .7 + nonstoich. spinel, ac, A1 2 C>3, C glass, Fe, be, graphite, hibonite, Fe30 4 

0.004 ± 0.001 

Ne II 

A1 2 0 3 + be, ac, C glass, Fe, graphite, FeS, Fe30 4 , Mgo.1Feo.9S 
(nonstoich. spinel + be, C glass, Fe) 

0.004 ± 0.001 

X-ray Fe 

Mg 2 . 4 Si0 4 .4 + be, C glass, Fe, Mgo.1Feo.9S, ac, graphite, FeS, Fe30 4 
MgFeSi0 4 + be, C glass, Fe, Mgo.1Feo.9S, ac 
(other silicates + usual “featureless” components) 

0.02 ±0.01 

South Spot 

MgFeSi0 4 + C glass, ac, be, Fe 

Mg 2 . 4 Si0 4 . 4 + ac, Fe, be, C glass 

(other silicates ± usual “featureless” components) 

0.0001 ±3 x 10 -5 

Radio 

Mg 2 . 4 Si0 4 . 4 ± graphite, nonstoich. spinel, Mgo.1Feo.9S, C glass, PAH+, Al 2 03, FeS, Fe, be 
MgFeSi0 4 + graphite, C glass, A1 2 C>3, Fe, be, FeS, Mgo.1Feo.9S 
(Mg 2 Si0 4 + C glass, ADO3, graphite, Fe, be) 

0.0004 ± 0.0004 

Si ii 

Indeterminate 

< 0.1? 


1 The compositions in bold need to be paired with any of the (generally) more featureless components listed after the “+” sign. These 
components are listed in order of increasing y 2 , but the differences are not significant. Combinations listed in parentheses have y 2 values 
larger than other combinations, but still within a factor of 2 of the minimum y 2 . 

2 Mean and a for good models in which the bold components at left are the dominant mass. 


Table 4. 


Possible Compositions for the 12 pm Peak of Ar II, Ar III Dust 


Compositions 

Comment 

Si0 2 

SiC 

Mg a ,Alj / 04 or CaAli 20 ig 

requires CDE calculation of also contributes to 9 and 21 /im features 

requires CDE calculation of k \ 

requires extremely high temperatures [y = (8 — 2x)/3] 


